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ABSTRACT 

We numerically investigate whether and how gaseous ejecta from AGB stars can 
be converted into new stars within originally massive star clusters (MSCs) in order to 
understand the origin of multiple stellar populations in globular clusters (GCs). We 
adopt a scenario in which (i) MSCs with masses of M s can be formed from high-mass, 
high-density giant molecular clouds (GMCs) in their host galactic building blocks 
embedded in dark matter halos at high redshifts and (ii) their evolution therefore can 
be significantly influenced by M s , their initial locations, and physical properties of 
their hosts. Our 3D hydrodynamical simulations show that gaseous ejecta from AGB 
stars can be retained within MSCs and consequently converted into new stars very 
efficiently in the central regions of MSCs, only if M s exceed a threshold mass (-Mth) 
of « 1O 6 M0. The new stars can correspond to the "second generation (SG)" of stars 
with higher Na and lower O abundances observed in GCs. Star formation efficiencies 
during the formation of SG stars within MSCs with M s ^ M tn can be rather high 
(0.3 — 0.9) so that very compact new clusters within original MSCs can be formed. M s 
should be as large as 10 6 — 1O 7 M to explain the observed large fraction of SG stars 
in the present ordinary Galactic GCs, because new stars can consist of only 1 — 4% 
among all stars for the standard IMF. Nuclear MSCs are found to retain much more 
effectively the AGB ejecta and convert more efficiently the gas into new stars owing 
to much deeper gravitational potential of their hosts. Capture and accretion of cold 
molecular gas (or small GMCs) by forming MSCs themselves can be mechanisms for 
mixing (i.e., dilution) of AGB ejecta with cold pristine gas. We suggest that both M s 
and their locations within their hosts can determine whether abundance spread can 
be seen only in light elements or even in heavy ones. We discuss how and in what 
time scale MSCs preferentially lose old stars owing to tidal stripping by their host 
galactic building blocks. We also suggest that the origin of the intermediate-age GCs 
with possible age spread of ~ 100 Myr yet apparently no/little abundance spread in 
light elements in the LMC is closely associated with their incapability to retain the 
AGB ejecta owing to their low masses. 
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1 INTRODUCTION 

It has long been discussed both observationally and theo- 
retically why some of the Galactic GCs show star-to-star 
inhomogeneity among the light elements of stars and what 
physical mechanisms are responsible for the inhomogeneity 
(e.g., Cottrell & Da Costa 1981; Sneden et al. 1992; Norris 
& Da Costa 1995; Cannon et al. 1998; Gratton et al. 2004; 
D Antona & Caloi 2004; Fenner et al. 2004; Norris 2004; Lee 
et al. 2005; Smith et al. 2005; Bekki et al. 2007; Alves-Brito 
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et al. 2008; Calelan 2008; Kayser et al. 2008; Piotto 2008; Da 
Costa et al. 2009; Marcolini et al. 2009; Yong et al. 2009; 
Carretta et al. 2010; D'Ercole et al. 2010; Romano et al. 
2010; van Loon 2010). Although such star-to-star inhomo- 
geneity was discovered in 70's and 80's (e.g., Cohen 1978; 
Peterson 1980; Norris et al. 1981; Leep et al. 1986), recent 
statistical studies for a larger number of the Galactic GCs 
have established that the presence of multiple stellar pop- 
ulations is an universal phenomena seen across most of the 
Galactic GCs (e.g., Carretta et al. 2010). The latest obser- 
vational results by Ferraro et al. (2009) have revealed that 
the Galactic metal-rich GC Tarzan 5 has two different stel- 
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lar populations with different abundances of heavy elements. 
Lee et al. (2009) has also suggested a significant fraction of 
the Galactic GCs have two different populations with differ- 
ent abundances in heavy elements based on color-magnitude 
diagrams of the GC stars in /ifc-bands. 

Now large star-to-star abundance variations have been 
confirmed in almost unevolved stars in some of the Galactic 
GCs (e.g., Gratton et al 200f ; Ramirez & Cohen 2002; Bedin 
et al. 2004; Gratton 2004; Carretta et al. 2004; Piotto et al. 
2005, 2007; D'Orazi et al. 2010a), which strongly suggests 
that after the first generation (FG) of stars formed within 
forming GCs, gas chemically mixed with gaseous ejecta of 
FG stars was converted into to the second generation (SG) 
of stars (and even third and fourth generations). Observa- 
tions have now being extensively investigating how chemical 
properties of SG stars (e.g. fractions of SG stars) correlate 
with their global internal properties (e.g., magnitudes and 
ellipticities) and with locations and 3D motions with respect 
to the Galactic center to understand the origin of SG stars 
and physical relationships between formation of FG and SG 
stars (e.g., Carretta 2006, Carretta et al. 2010). 

Following these observational developments, theoretical 
studies have considered that SG stars can form from gaseous 
ejecta either by AGB stars ("AGB scenario"; D'Antona & 
Caloi 2004; Karakas et al. 2006) or by fast rotating massive 
stars ( "FRMS" scenario; Prantzos & Charbonnel 2006; De- 
cressin et al. 2007) and thereby investigated whether the 
observed physical properties of the Galactic GCs can be 
explained by their models. Although it has not been de- 
termined which of the two can explain observations in a 
self-consistent manner, a number of authors suggested that 
the AGB scenario is a more physically viable because AGB 
ejecta is more likely to be converted into new stars owing to 
low ejection velocities (e.g., Renzini 2008). 

Recent observations have provided the following two 
key results which can give strong constraints on any the- 
ory for GC formation: (i) the large fraction (typically 0.67) 
of SG stars in each individual Galactic GCs and (ii) the 
observed Na-0 and Mg-Al anti-correlations between cluster 
stars (e.g., Carreta et al. 2010 for a comprehensive study 
for these two). The first key result suggests that the orig- 
inal stellar systems (i.e., FG stars) need to be much more 
massive than the present GCs unless we adopt very unusual 
and unrealistic IMFs (e.g., Smith & Norris 1982; DAntona 
& Caloi 2004): the original systems are either very massive 
star clusters or dwarf galaxies hosting GCs (e.g., Bekki & 
Norris 2006). If the original systems are really massive ones, 
then the FG stars need to be preferentially lost while SG 
ones remain the same so that the observed typical fraction 
(~ 0.67) of SG stars can be explained. 

Although a number of theoretical works based on the 
AGB scenario tried to explain the second key observational 
result in a self-consistent manner (e.g., Fenner et al. 2004; 
Bekki et al. 2007; D'Antona et al. 2005; D'Antona & Ventura 
2007; Ventura & D'Antona 2006, 2008; D'Ercole et al. 2010), 
their models appear to have not yet explained all of the rel- 
evant observations on chemical abundances of the Galactic 
GCs in a fully self-consistent manner. Chemical evolution 
models based on the FRMS scenario have not yet been fully 
explored so that the validity of the FRMS scenario can not 
be currently assessed. 

In these previous models, it is assumed that AGB ejecta 



can be converted into new stars (i.e., formation of SG stars) 
within already existing clusters. However it would not be 
so obvious that such secondary star formation can occur 
within clusters, given the shallow gravitational potential 
wells of clusters and the small mass fraction of AGB ejecta. 
Therefore, secondary star formation processes within clus- 
ters should be investigated by numerical simulations that 
can include various physical processes within clusters (e.g., 
retention of AGB ejecta). D'Ercole et al. (2008) first inves- 
tigated whether SG stars can be formed from the gaseous 
ejecta of AGB stars of FG with M B = 10 7 M Q . However, their 
models are based on one-dimensional hydrodynamical simu- 
lations and have limitations in predicting 3D structures and 
kinematics of final stellar systems. Three-dimensional (3D) 
stellar and gas dynamical numerical simulations with a plau- 
sible model for star formation are ideal to investigate sec- 
ondary star formation within clusters and can furthermore 
provide theoretical predictions that can be compared with 
the observed differences in 3D structures and kinematics be- 
tween FG and SG stars (e.g., Norris et al. 1997; Ferraro et 
al. 2002; Sollima et al. 2005, 2007; Pancino et al. 2007 Bellini 
et al. 2009; Anderson & van der Marel 2010). 

The purpose of this paper is thus to investigate exten- 
sively star formation from gaseous ejecta from AGB stars 
within MSCs based on self-consistent hydrodynamical simu- 
lations with a reasonable model for star formation. We adopt 
a scenario in which (i) MSCs can be formed in their host 
galactic building blocks embedded in dark matter halos at 
high redshifts and (ii) new stars formed from AGB ejecta 
can finally become the observed SG stars in the present 
GCs. Based on this scenario, we investigate (i) how AGB 
ejecta can be retained within MSCs, (ii) whether and how 
secondary star formation from AGB ejecta proceeds within 
MSCs, and (iii) how MSCs evolve if they are located in nuclei 
of their hosts. We also investigate how MSCs lose their stars 
during tidal interaction with their hosts in order to discuss 
the observed smaller fractions of FG stars in the Galactic 
GCs. 

A number of previous works discussed a scenario in 
which GCs were originally stellar galactic nuclei or nuclear 
star clusters in nucleated galaxies and their host galaxies 
had been already destroyed by strong tidal fields of much 
larger galaxies to disappear completely (e.g., Zinnecker et 
al. 1988; Freeman 1993; Bekki & Freeman 2003; Bellazzini 
et al. 2008; Boker 2008). Following this scenario, Bekki et 
al. (2007) investigated chemical abundances of FG and SG 
stars in GCs formed within the central regions of their host 
galaxies. The evolution of MSCs initially nuclear regions of 
their hosts in the present study therefore can provide some 
implications on the validity of the above scenario in explain- 
ing observational properties of GCs. 

The plan of the paper is as follows: In the next sec- 
tion, we describe the details of the proposed scenario for 
GC formation. In §3, we describe the numerical models for 
evolution of gaseous ejecta of AGB stars within MSCs. In §4, 
we present the numerical results mainly on physical proper- 
ties of SG stars. In §5, we discuss how long MSCs can retain 
most of their stars when they are influenced by strong tidal 
fields of their hosts. In §6, we discuss a number of key issues 
related to the origin of multiple stellar populations in GCs. 
We summarize our conclusions in §7. 
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Figure 1. Dependences of the ratio (/ e j) of the total mass of 
AGB ejecta (JWagb) to the total stellar mass of a MSC (M s ) 
as a function of the lower- mass cut-off of AGB stars (mi agb) 
and the time that has elapsed since the most massive AGB star 
(rrt u agb = 8Mq) starts to eject gas for two different models with 
the IMF slopes of a = 2.35 (blue solid) and 1.5 (red dotted). This 
/ c j describes a possible maximum mass fraction of AGB ejecta 
that can be converted into new stars within a MSC. Note that 
at most ~ 10% of Ms can be possibly converted into new stars 
within 1 Gyr. 

2 THE SCENARIO 

Although we describe the scenario in the context of the 
Galaxy formation, the formation processes of GCs within 
galaxies in general would be similar to those described be- 
low. The scenario is based on the results from our present 
and previous theoretical studies on GC formation (e.g., 
Bekki et al. 2002; Bekki et al. 2004; Bekki & Chiba 2007; 
Bekki 2006; Bekki et al. 2007; Hurley & Bekki 2008; Bekki et 
al. 2008). Since the scenario is based partly on results of cos- 
mological simulations (e.g., Bekki et al. 2008), the scenario 
has some implications on the observed correlations between 
abundance properties of GCs and 3D motion and kinematics 
of GCs (e.g., Carretta 2006; Lee et al. 2007). In the present 
paper, the masses of the present GCs are represented by 
M gc so that M s (original masses of the GCs) and Af gc can 
be discriminated with each other. 

2.1 Time sequence 

2.1.1 Formation of FG stars of MSCs within host dwarfs 

FG stars of MSCs form from high-density GMCs with 
masses (Mgmc) as large as or larger than 1O 7 M0 within 
the Galactic building blocks which are later destroyed ow- 



ing to the strong tidal field of the Galaxy when they merge 
with the Galaxy. Most of the hosts form as massive dwarf 
galaxies with masses (Mh) larger than 1O 9 M0 before reion- 
ization (z > 10), and thus their star formation and chem- 
ical enrichment histories can be different from those of the 
present dwarfs within the Galaxy (Bekki et al. 2008). The 
formation places of MSCs are highly likely to their hosts' 
central regions where mass fractions of GMCs can be much 
higher (i.e, not necessarily in the nuclei of the hosts). 

The MSCs formed in the very centers of their hosts 
start their lives as nuclear MSCs so that their evolution can 
be different significantly from that for MSCs formed outside 
the very centers owing to much deeper gravitational poten- 
tials of hosts' central regions. The host GMCs of MSCs may 
well initially have numerous substructures (smaller GMCs) 
owing to their large masses (e.g., Efremov 1995; Bonatto & 
Bica 2010) so that MSCs at their birth can be composed of 
a number of smaller clusters (i.e., not single entities). These 
smaller sub-units can finally merge with one other to form 
single massive clusters within the merging timescale depend- 
ing on M s and their sizes (Rs). 

2.1.2 Evolution of FG stars within GMCs 

Massive stars and type II supernova can strongly influence 
later evolution of molecular gas left behind from the forma- 
tion of FG stars. If most of the energy from massive stars 
and type II supernova are converted into kinetic energy of 
gas surrounding the stars, then gaseous ejecta from the stars 
can not be mixed well with the residual molecular gas (Bekki 
& Chiba 2007). Therefore, SG formation from gaseous ejecta 
from the above energetic stars and the mass fraction of SG 
stars formed from the mixed gas is very small (< 10 -4 ). 
However, if massive stars are FRMSs with very slow ejec- 
tion radial velocities (< 10 km s _1 ), then the gaseous ejecta 
could be well mixed with the residual gas to form SG stars 
before supernova explosion expels a significant fraction of 
gas. 

Although type II supernova can heat up and expel a sig- 
nificant fraction of gas from massive GMCs, some fraction 
of gas can be immune from the effects of supernova owing 
to the clumpy nature of the massive GMCs. The smaller 
residual clouds may well be later accreted onto the forming 
clusters to dilute AGB ejecta. While supernova explosion is 
influencing significantly GMCs, star formation rates within 
the GMCs can drop significantly. Gaseous ejecta from mas- 
sive AGB stars start to be accreted onto the central regions 
of MSCs after all supernova events finish. However star for- 
mation does not resume until an enough amount of gas can 
be accumulated within the central regions: there should be a 
threshold gas mass fraction (/ e j,th) above which star forma- 
tion can start within MSCs. Therefore there can be signif- 
icant time-delay (~ 10 s yr) between FG and SG formation 
in this scenario. 



2.1.3 Formation of SG stars 

SG stars in a MSC start to form efficiently in the central 
region of the MSC when the mass fraction of the accumu- 
lated AGB ejecta to M s (f e j) exceeds / e j,th- This efficient sec- 
ondary star formation from AGB ejecta can occur only if M s 
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Table 1. Range of model parameters for numerical simulations of star-forming MSCs 

M a ^_ R s £ StotJ; Iagb£ Star formation^ p th J_ External tidal field 3_ M h £ 

10 4 -To 8 35 - 200 - 0~2 0.002 - 0*08 YES/NO ' - 10 4 YES/NO 10 8 - 2 x~10 1( 

a The initial mass of a MSC in units of Mq. 
^ The initial size of a MSC in units of pc. 

c The ratio of rotational energy to total kinetic one in a MSC. 

^ The ratio of the total mass of AGB stars to that of stars in a MSC. 

e "YES" ("NO") means that the star formation model is (is not) included in the simulation. 
f The threshold gas density for star formation in units of atoms cm -3 

9 "YES" ("NO") means that a MSC can (can not) be influenced by the tidal field of its host galaxy. 
^ The total mass of the host galaxy for a MSC in units of Mq. 



exceeds a threshold cluster mass (M t h), because AGB ejecta 
can not be efficiently retained in MSCs with smaller M s . 
During this secondary star formation, the residual molecu- 
lar gas can be mixed with AGB ejecta to form new stars. 
Owing to the presence of deeper gravitational potential wells 
of MSCs, star formation efficiencies (e s f ) can be much higher 
(> 0.3) for this secondary star formation so that very com- 
pact clusters can be formed. MSCs thus can initially show 
"nested structures" with more diffuse distributions of FG 
stars and more compact ones of SG stars. 

Although MSCs with M s below M t h can hardly form 
new stars from AGB ejecta, such less massive clusters can 
still capture residual molecular gas and smaller nearby 
GMCs (not chemically contaminated by ejecta from FG 
stars). Cold molecular gas with initially very high-densities 
obtained by MSCs themselves can be converted into new 
stars within the central regions of MSCs (Bekki & Mackey 
2009). Very low- mass MSCs can be destroyed by interaction 
with GMCs so that secondary star formation can not hap- 
pen in their central regions. Therefore, there should be a 
threshold cluster mass below which MSCs show differences 
neither in ages nor in light elements among their stellar pop- 
ulations. Secondary star formation processes within nuclear 
MSCs can be significantly different from those described 
above owing to much deeper potential wells of their hosts 
(as described in the present study). 



2.1.4 Tidal stripping of MSCs due to their host dwarfs 

Secondary star formation can continue until MSCs lose most 
of their FG stars, though massive stars and type II super- 
nova of SG stars can also suppress or even truncate further 
star formation. Strong tidal fields of hosts can efficiently 
strip stars preferentially from FG stars owing to initially 
diffuse spatial distributions. This preferential stripping by 
the Galaxy has been already proposed by D'Ercole et al. 
(2008) and discussed in the context of origin of the Galac- 
tic stellar halo (Vesperini et al. 2010). The present scenario 
suggests that not the tidal field of the Galaxy but those of 
hosts of MSCs are responsible for the tidal stripping of most 
FG stars in MSCs. The timescale of a MSC to lose most of 
their FG stars due to tidal stripping by its host depends 
on its position, with respect to its host (R p ), M B , and R B . 
Therefore, more massive and denser MSCs can continue sec- 
ondary star formation longer so that AGB stars with lower 
masses can possibly contribute to further star formation. 



2.1.5 Evolution into the Galactic halo GCs 

Hosts of MSCs are strongly influenced by the tidal field of 
the forming Galaxy during their merging with the Galaxy 
so that they can be completely destroyed by the Galaxy. 
The stripped stars from hosts, which include FG stars from 
MSCs, can be some parts of the Galactic stellar halo. MSCs 
are stripped during disintegration of their hosts to finally be- 
come halo GCs. Strong ram pressure of the Galactic halo gas 
(e.g., Frank & Gisler 1976; Bekki 2006) and no cold gas avail- 
able in the halo prevent almost completely star formation of 
the MSCs since they become the Galactic halo GCs. Thus, 
multiple stellar populations can be formed within MSCs only 
when they are within their hosts. 

Gas can be efficiently transferred to the nuclear regions 
of hosts to be converted into new stars while they are be- 
ing destroyed by the Galaxy (e.g., Bekki & Freeman 2003). 
These new stars can have chemical abundances in heavy ele- 
ments different from those of original nuclear MSCs, because 
they are from the outer regions, where chemical enrichment 
histories are quite different from those of nuclei. Nuclear 
MSCs can show star-to-star abundance variations not only 
in light elements but also in heavy ones. 

2.2 Constraints on original masses of MSCs 

A number of authors have already suggested that if the stan- 
dard IMF is applied for FG stars and if e s f (star forma- 
tion efficiency for SG formation) is 1.0, then original total 
masses of FG stars (M s ) are required to be at least ten times 
larger than the present masses of SG stars in the Galactic 
GCs with multiple stellar populations (e.g., Bekki & Nor- 
ris 2006). Given that e s f is not like 100% as observed in 
star-forming regions in the Galaxy, the required M s should 
be even higher than the above. Recent chemical evolution 
models have shown that the total mass of pristine gas that 
can mix with AGB ejecta to form SG stars needs to be 
comparable to that of AGB ejecta in order to explain the 
observed levels of star-to-star abundance variations and the 
O-Na anti-correlation in the Galactic GCs (e.g., D'Ercole 
et al. 2010): AGB ejecta should not be too much diluted 
by external gas captured/accreted by MSCs. Thus the re- 
quired M s estimated in previous studies can not change sig- 
nificantly even if contribution from external pristine gas to 
SG formation is considered. 

In the present scenario, gaseous ejecta from AGB stars 
with masses ranging from tti^agb to m u ,AGB can contribute 
to the formation of SG stars within a timescale of t s (which 
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determines iti^agb)- We here estimate the mass fraction of 
AGB ejecta to the initial mass of a MSC (/ e j) (i) for a given 
mi,AGB and (ii) for a given t s by assuming m u ,AGB = 8M© 
and power-law initial mass functions with the slopes of a. 
The adopted IMF in number defined as ip(mi) = M St omi~ a , 
where mi is the initial mass of each individual star and the 
slope a = 2.35 corresponds to the Salpeter IMF. The nor- 
malization factor M St o is a function of M s , m\ (lower mass 
cut-off), and m u (upper mass cut-off): 

_ M s x(2~a) 

where mi and m u are set to be O.1M and 12OM , respec- 
tively. The total mass of AGB ejecta within a MSC (Magb) 
is accordingly described as: 

/■™ii,AGB 

Magb = / m ei ip(m)dm, (2) 

m l,AGB 



T = (Myr) 



2.7 



where m c j describes the total gas mass ejected from an AGB 
star with initial mass mi and final mass (mp). We derive an 
analytic form of m e j (= mi — mp) from the observational 
data by Weidemann (2000) by using the least-square fitting 
method, and find: 



m ei = 0.916Mi - 0.444. 



(3) 



In order to calculate the main-sequence turn-off mass 
(triTo) we use the following formula (Renzini & Buzzoni 
1986): 



logm T o(t s ) = 0.0558(logi s ) - 1.338 log t s + 7.764, 



(4) 



where mxo is in solar units and time t s in years. 

Fig. 1 shows that / c j is larger for smaller mi,AGB both 
for the standard IMF (a = 2.35) and the top-heavy one 
(q = 1.5). Although the top-heavy model shows larger / e j 
than the standard one for mi,AGB < 2M , the derived / c j 
is small: only 0.036 for rn^AGB = 5M and 0.079 rn^AGB = 
3M . Ventura & D'Antona (2008) showed that the observed 
O-Na anti-correlations can be well reproduced only if AGB 
stars with masses larger 5M can contribute to secondary 
star formation within GCs. Thus, it is highly likely that 
progenitor clusters for the present GCs. needs to be at least 
~ 25 times more massive than the present GCs. 

As shown in Fig. 1, / e j can be larger than 0.1 for t s > 1 
Gyr (or mi,AGB < 3M ) for the two IMF models. It should 
be here stressed that MSCs with larger t s can end up with 
a large age difference (> 1 Gyr) between FG and SG stars. 
Original MSCs can lose most of their FG stars due to tidal 
stripping of their hosts within well less than ~ 1 Gyr, as 
shown later in the present study. Therefore gaseous ejecta 
only from massive AGB stars can participate secondary star 
formation within MSCs. It would be reasonable to consider 
that mi j agb is 4 — 5M that corresponds to t s ~ 10 8 yr. 



3 THE NUMERICAL MODEL 

We investigate star formation of gaseous ejecta from FG 
AGB stars in MSCs with M s ranging from 1O 4 M to 1O 8 M 
by using the latest version of GRAPE (GRavity PipE, 
GRAPE- 7) which is the special-purpose computer for grav- 
itational dynamics (Sugimoto et al. 1990). We have revised 
our original GRAPE-SPH code (Bekki 2009) so that we 
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Figure 2. Time evolution of the distribution of gaseous ejecta 
from an AGB star with V e j = 20 km s _1 projected onto the x-z 
plane for models with M s = 2 X 1O 5 M (upper four panels) and 
M s = 1O 6 M (lower four). 



can investigate star formation processes within the above- 
mentioned massive stellar system. MSCs with very large 
(M s > 1O 7 M ) are still referred to as "clusters" just for 
convenience in the present study, though their stellar masses 
are as massive as those of dwarf galaxies 

We focus mainly on secondary star formation from 
AGB ejecta of FG stars within original MSCs and describe 
the results of numerical simulations on the star formation. 
We have numerically investigated ram pressure stripping of 
AGB ejecta by ISM of their hosts and its strength dependent 
on M s and physical parameters of ISM. We however describe 
the results in our forthcoming papers (Bekki et al. 2010), 
because ram pressure stripping of AGB ejecta of MSCs by 
typical ISM of hosts is only effective for low-mass MSCs with 
M s < 1O 4 M , (which can not become ordinary GCs in the 
present study) and thus not so important. 



6 K. Bekki 



8xl0" 5 



a 6xl0" s 



4xl0" 6 



E 2xl0" s 



-i— i— i— i— p-i— I— r 

F„ t (R<50pc) 
F rel (R<10pc) 



4 4.5 5 5.5 6 

log 10 M a (M ) 

Figure 3. Dependences of the mass ratios (.Fret = Magb/Ms) 
of AGB ejecta that can finally be within the central 50pc (blue 
filled circles) and lOpc (red filled squares) of MSCs on M B . 



We have numerically investigated accretion of warm 
ISM onto MSCs by Bondi accretion and found that Bondi 
accretion with the accretion rate of ~ 1O _3 M yr" 1 is possi- 
ble only if M s is larger than 3 x 1O 6 M for typical ISM with 
n — 1 atom cm -3 , relative velocity of ~ 20 km s~\ and gas 
temperature of 10000 K. Although this accretion rate is too 
small for MSCs to obtain fresh gas for the formation of SG 
stars within ~ 10 s yr, this Bondi accretion can be an im- 
portant mechanism for obtaining fresh gas in massive GCs 
with abundance variations in heavy elements like cj Cen, as 
briefly discussed later in this paper. We will describe the de- 
tails of the numerical results on Bondi accretion of ISM in 
our forthcoming papers (Bekki et al. 2010), mainly because 
this paper becomes too long if the results are included in 
this paper. 



3.1 Initial structures and kinematics of MSCs 

The original MSCs are assumed to have a Plummer den- 
sity profile (e.g., Binney & Tremaine 1987) with luminosities 
(L sc ) and central velocity dispersions (er s ) consistent with 
the relation observed for GCs (Djorgovski et al. 1997): 



K a s 



(5) 



where Ko is a normalization factor for the relation. The scale 
length (a s ) of a MSC is determined by the formula 



a B — GM B /6a, 



2 
s ) 



(6) 



where G is the gravitational constant. Since the-mass-to- 
light-ratio (M sc /L s ) is assumed to be constant for all SCs, 
a s and <r a are determined by the equations (1) and (2) for 
a given M sc . The normalization factor Ko in equation (5) 
is determined such that a cluster with M s — 6 x 10 Mq 
can have the size of R s = 50pc (R s = 5a s ), and the central 
velocity dispersion of 7 km s" 1 . 

Hasegan et al. (2005) revealed that massive stellar sys- 
tems with masses larger than 2 x 10 6 Mq appear to have 
scaling regions different from those of the present GCs. 
We accordingly consider that original sizes of MSCs with 
M s ^ 2 x 1O 6 M can not be determined by the above equa- 
tions (5) and (6) and therefore investigate models with dif- 
ferent Rs for a given M s for such MSCs. It is also possible 
that original MSCs with M s < 2 x 1O 6 M have 7? s smaller 
or larger than those described by the above scaling relations 



because MSCs composed only of FG stars have ages of ~ 10 
yr in the present study: an order of 10 8 yr would not be long 
enough to form GCs on the present scaling-relation owing 
to dynamical relaxation processes. We therefore investigate 
models with i? s smaller and larger than the above equation 

(5) and (6) predict for a given M s . 

A MSC is assumed to have a small amount of rotation 
with the ratio of the initial rotational energy (Trot) to the 
total kinetic one (T^ n ) being a free parameter represented 
by Srot- The parameter values of s rot range from (no rota- 
tion) to 0.3 (rapid rotation). The initial rotational velocity 
of a particle at a distance of R from the center of MSC is 
cj R, where ui (constant angular velocity) is determined such 
that s r ot can be the adopted value. Therefore, the system 
has random kinetic energy (T ra n) of (1 — s ro t)2ki n (due to 
isotropic velocity dispersion of stars). Firstly we estimate 
Tkin for a s determined by M s and Rs (in equations (5) and 

(6) ) and then reduce u s so that the final system can be in 
virial equilibrium (Tkm = T ran + Trot = 0.5|W|, where W is 
the total potential energy of the system). 

3.2 Gas dynamics 

The FG stars in a MSC are represented by equal-mass stellar 
particles with the particle number of N s (= 10 ) and some 
fraction of the particles are assigned as "AGB stars" with 
initial masses of m s and the total number of A?agb that can 
eject SPH gas particles with ejection velocities of V e j with 
respect to the centers of the AGB stars. The mass fraction 
of AGB particles in a MSC is a parameter represented by 
/agb that is determined by the adopted IMF. As outlined 
in §2, gaseous ejecta only from massive AGB stars with ini- 
tial masses of ~ 5 — 8M need to be converted into new 
stars (i.e., SG stars) to explain the chemical abundances of 
SG (e.g., Ventura & D'Antona 2008). In the present paper, 
V e j=20 km s is adopted, which is a reasonable choice for 
massive AGB stars (e.g., Marshall et al. 2004). 

The present simulations can not resolve gaseous evo- 
lution of each individual AGB star owing to the adopted 
numerical resolution (~ 0.5 pc). We therefore assume that 
each AGB particle initially has an expanding gaseous sphere 
which is much larger than the AGB star itself. The mass, 
size, and temperature of the large gaseous sphere repre- 
sented by SPH particles with the particle number of tiagb 
are set to be m g , r g and T g , respectively. Each AGB parti- 
cle accordingly has a gas sphere represented by tiagb SPH 
particles with radial velocities (with respect to the AGB 
particle) of V j. The AGB particle therefore has a mass of 
m s — m g after gas ejection. We here assume that each AGB 
particle in a model eject SPH particles at T = when the 
simulation starts: we do not consider gradual ejection of gas 
from AGB stars at each time step, because an appropriate 
modeling of such gradual ejection requires a huge number of 
gas particles and thus is numerically costly (i.e., practically 
not feasible). 



Although we have investigated models with T g 



100- 



1000K, we mainly describe the results of the models with 
T g = 100K which corresponding to star-forming warm 
molecular clouds (e.g., Wilson et al. 1997). This is because 
we consider that AGB wind can cool down during expan- 
sion throughout interstellar space due to radiative cooling 
(to finally become molecular gas for further star formation) 
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4. The same as Fig. 2 but for the model with two AGB 
a MSC with M s = 2 x 1O 5 M . 



so that T g can becomes much smaller than the original tem- 
perature of the wind (~ 1000K). We adopt an isothermal 
equation so that AGB ejecta can keep initially low temper- 
atures (100K). It should be stressed here that if T = 1000K 
is adopted (which is not realistic though), star formation is 
possible only for models with M s ^ 10 6 Mq. 

Most of AGB ejecta can be converted into new stars 
well before a few Myr after commencement of secondary star 
formation so that feedback effects of massive stars and type 
II supernova can not significantly influence gas dynamics 
during the secondary star formation processes. We therefore 
consider that the above isothermal assumption is reasonable. 
In order to estimate value of m g corresponding to the total 
mass ejected from each AGB star after the main-sequence 
turn-off, we use the formula given in the equation (3). About 
83% of a AGB particle with m s = 5M© can be ejected to be 
used for further star formation in the present model. 

If we adopt the Salpeter IMF, then /agb (the mass 
fraction of AGB stars with masses larger than 5Mq) is ~ 
0.04 in the present study. Since initial masses within a MSC 
are the same between stellar particles, Nagb = /agb-W s . 
Although we adopt this /agb = 0.04 for most models, we 
investigate models with different /agb so that we can find a 
threshold /agb for a given MSC above which secondary star 
formation can occur in the MSC. It is found that if /agb 
is larger than 0.008, secondary star formation is possible 
in some models, though the star formation efficiency is low. 
This means that secondary star formation within a MSC can 
not be possible until a certain amount of gas is accumulated 
within the MSC and thus that there should be a time-delay 
between the formation of FG stars and that of SG ones. 



3.3 Star formation 

We investigate whether the gas accumulated in the central 
regions of MSCs can be sufficient to form new stars by adopt- 
ing a simple prescription for star formation. In the mod- 
els with "star formation" , a gas particle is converted into a 
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collision-less new stellar one if the gas particle meets the fol- 
lowing three conditions: (i) the dynamical time scale of the 
SPH gas particle is shorter than the sound crossing time, 
(ii) the gas is converging (i.e., Vv < 0, where v is the veloc- 
ity vector of the gas particle), and (iii) the local gas density 
exceeds a threshold gas density (pth) which corresponds to 
the densities of dense cores of molecular clouds (n > 10 4 
atoms cm -3 ), where individual star formation is ongoing. 
The first two conditions mimic the Jeans gravitational in- 
stability for gaseous collapse. The new stellar particles and 
old ones initially within MSCs are referred to as "new stars" 
(i.e., SG stars) and "old stars" (i.e., FG ones), respectively. 
Also original clusters and new ones formed from AGB ejecta 
are referred to as "old clusters" and "new clusters" , respec- 
tively, just for convenience. 

We investigate models with p t h = and 10 4 atoms cm -3 
in order to investigate (i) whether the threshold gas den- 
sity is important for determining star formation histories 
of SG stars and (ii) how the present results, in particular, 
structural and kinematical properties of final clusters com- 
posed of SG stars, depend on p t h: we here consider that pth 
has not been observationally well determined yet and there- 
fore could be different for individual star forming clouds. 
Although these star formation models are less realistic in 
some points (not inclusive of magnetic fields and radiative 
transfer etc), we consider that they enable us to grasp es- 
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Figure 7. Time evolution of the distributions of gas (magenta) and new stars (cyan) projected onto the x-z plane for the standard 
model with M s = 10 6 M©, R s = 44.7pc, s r ot = 0.0, and p t h = 10 4 atoms cm -3 . The green circle in each frame represent the original 
half-mass radius of old stars. Time (T in units of Myr) that has elapsed since the simulation started is shown in the upper left corner 
for each frame. 



sential ingredients of secondary star formation from AGB 
ejecta within original massive clusters. 



3.4 External gravitational fields of host galaxies 

We mostly investigate "isolated MSCs" which are located 
outside nuclei of their host galaxies and therefore short-term 
(~ 10 7 yr) evolution of AGB eject can not be strongly influ- 
enced by the tidal fields of their hosts. However, dynamical 
evolution of AGB ejecta from nuclear MSCs that are initially 
located in the very centers of their hosts can be significantly 
influenced by their hosts owing to the much deeper gravita- 
tional potential wells of the hosts. We therefore investigate 
how hosts influence evolution of AGB ejecta from nuclear 
MSCs by assuming that the host dwarfs are dominated by 
dark matter halos with masses of Mh and thus their gravi- 
tational potentials are determined only by the dark halos. 

We adopt the density distribution of the NFW halo 
(Navarro, Frenk & White 1996) suggested from CDM simu- 
lations: 

p(r) = £2 (7) 

P[> (r/r s )(l + r/r s ) 2 ' (> 

where r, po, and r a are the spherical radius, the characteristic 



density of a dark halo, and the scale length of the halo, 
respectively. Although we investigate models with 10 8 Mq $J 
Mh ^ 2 x 1O 1O M0, we mainly describe the results of the 
models with Mh = 10 9 Mq, r s = 0.6 kpc, and the virial 
radius (r v - a ) of 9.9 kpc (i.e., the c-parameter of 13.9). These 
models with Mh = 1O 9 M0 are reasonable in the present 
scenario in which there is a threshold galaxy mass (1O 9 M ) 
above which GCs can be formed. It should be here stressed 
that low-mass dark matter halos with Mh = 1O 8 M0 can 
influence evolution of AGB ejecta in hosts as significantly 
as those with M h = 1O 9 M . 



3.5 Ranges of parameters 

Firstly we investigate (i) whether AGB ejecta can be re- 
tained within MSCs and (ii) whether and how secondary 
star formation from AGB ejecta can proceed for models that 
are consistent with the scaling relation of GCs shown in the 
equation (2). Secondly we investigate structure and kinemat- 
ics of the simulated clusters in the models with M B — 1O 7 M0 
and R s = lOOpc in which the final masses of SG stars can be 
as large as ~ 10 5 Mq. Thirdly we investigate whether much 
deeper gravitational potential wells of host dwarf galaxies 
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can play a role in retaining AGB ejecta efficiently in the 
nuclear MSCs. 

We have investigated secondary star formation pro- 
cesses for numerous models with different model parame- 
ters (e.g., M a , _R S , p t h, and s rot ); we mainly describe here 
the results for the "standard model" with M s = 1O 6 M0, 
R s — 44.7pc, s rot = 0.0, and p tri = 10 4 atoms cm -3 . We as- 
sume that 7i agb = 20 for most models except those for inves- 
tigating how AGB ejecta can be retained within MSCs. The 
reason for this adoption is that as long as wagb > 10, the 
results do not depend on tiagb- Since we focus on secondary 
star formation processes with a timescale of < 10 7 yr, we do 
not intend to discuss long-term dynamical evolution of clus- 
ters due to two-body relaxation. We accordingly introduce 
a gravitational softening length (e g ) for each simulation and 
e g is set to be equal to the mean particle separation at the 
half-number radius of old stars. 

The ranges of model parameters are shown in the Table 
1. We describe in detail only the results of some representa- 
tive models showing key parameter-dependences in retention 
processes of AGB ejecta and secondary star formation pro- 
cesses within clusters among the investigated models. Sec- 
ondary star formation processes can continue only for ~ 10 7 
yr so that they can not be strongly influenced by gravita- 
tional fields of their hosts unless MSCs are located in the 
very centers of their hosts. We thus show mainly the results 
for models without gravitational fields of hosts in order to 
much more clearly show the importance of the four key pa- 
rameters M a , Rs, s rot , and p t h in secondary star formation 
processes. 



3.6 Limiations of the model 

As described in §3.2, we do not consider that different AGB 
stars with different masses eject gas at different T owing to 
different lifetimes of the stars. Therefore all of the gaseous 
ejecta can be converted into new stars within ~ 10 Myr. 
This rapid consumption of the AGB ejecta is not so realis- 
tic, given that there should be at least several tens millions 
years time delay between the epochs when AGB stars with 
masses of 8M© and 5Mq start to eject their gas. If we con- 
sider the different epochs of gas ejection from AGB stars 
with different masses, then the star formation period would 
become significantly longer. However, energy feedback ef- 
fects of numerous supernovae from the SG stars are likely to 
truncate star formation in gaseous ejecta from lower-mass 
AGB stars. In our future papers, we will discuss this point 
in detail by using a more sophisticated numerical model. 

The present model does not include the effect of type la 
supernova (SN la) on star formation processes of SG stars. 
D'Ercole et al. (2008) already showed that the cumulative 
effect of numerous SN la is important for the star formation 
processes within massive clusters, because it can drastically 
alter gas dynamics within the clusters. Therefore the present 
model that does not include such SN la effect possibly over- 
estimates the star formation rates in the accumulated AGB 
ejecta in the central regions of MSCs. Our future more so- 
phisticated models with SN la effect will discuss how SN la 
can influence gas dynamics and star formation in MSCs. 
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standard model. 
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parameters being the same as those of the standard model (red 
dotted). 
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4 RESULTS 

4.1 Retention of AGB ejecta 

Fig. 2 show that time evolution of spatial distributions of 
AGB ejecta is very different between models with M s — 2 x 
10 5 Mq and M B = 1O 6 M0 in which model parameters except 
Ms are exactly the same. Gaseous particles with jiagb = 10 4 
ejected from an AGB star locating in the very center of a 
MSC can expand quickly (T = 2.7Myr) and finally escape 
from the cluster and never return back in the model with 
M s = 2x 10 5 M Q (T = 13.6Myr). On the other hand, gaseous 
ejecta from an AGB star in the model with M s = 1O 6 M0 
can return back to the initial location within a timescale of a 
few Myr after the ejection owing to the deeper gravitational 
potential for this more massive cluster. This clear difference 
suggests that M s is an important parameter for determining 
whether AGB ejecta can be retained in MSCs during their 
early formation phases. 

Fig. 3 shows how the mass fractions of AGB ejecta re- 
tained within MSCs (= F Ict ) depend on M s . Fig. 3 shows 
that MSCs with M s > 2 x 1O 5 M can retain a significant 
fraction of AGB ejecta if a scaling-relation similar to the 
observed one is applied for all MSCs. This result implies 
that there can be a threshold mass above which AGB ejecta 
can be efficiently retained. For MSCs with M s ~2x 10 5 M Q , 
AGB ejecta initially in the outer part of the MSC can escape 
from the MSC whereas that initially in the inner part can be 
retained within the MSC. As a result of this, a "core-halo" 
structure of AGB ejecta can be formed within 13.6 Myr. A 
very compact gas sphere can be formed where star forma- 
tion would be able to proceed rapidly if star formation is 
included in the model with M s — 1O 6 M0. A larger fraction 
of AGB ejecta can be retained in MSCs with larger M s . 

Gaseous spheres composed of AGB ejecta can inter- 
act with one another within a MSC when numerous stars 
enter into AGB phase almost simultaneously. The present 
3D numerical simulations enable us to investigate whether 
and how hydrodynamical interaction between the gaseous 
spheres increase or decrease the total gas mass accumulated 
within the central region of the MSC. Fig. 4 shows that if 
two AGB stars with each having wagb = 10 4 are included 
in the model with M s = 2 x 1O 5 M0 (where no gas accumula- 
tion is seen in Fig. 2), the gas can lose kinetic energy owing 
to energy dissipation during hydrodynamical interaction of 
the gaseous spheres so that a non-negligible amount of gas 
can be accumulated within the cluster. About 3.6% of the 
initial AGB ejecta can be accumulated within the central 
50pc of the cluster within 13.6 Myr. 

Fig. 5 shows how F let /NAGB, where Aagb are the to- 
tal numbers of AGB stars, depends on TVagb: this normal- 
ized F Ict can measure the efficiency of retaining AGB ejecta 
within MSCs. It is clear from this figure that AGB ejecta can 
be more efficiently retained within clusters for models with 
larger Aagb. This is mainly because as the larger number 
of gaseous spheres (of AGB stars) hydrodynamically inter- 
act with one another, the larger amount of the gas can lose 
their kinetic energy to be accreted onto the central regions 
of MSCs. The results in Figs. 4 and 5 therefore demonstrate 
that hydrodynamical interaction of AGB ejecta can cause 
efficient gaseous dissipation and thus play a great role in 
retaining AGB ejecta within MSCs. 

Fig. 6 shows that (i) there is a threshold M s (= 2 x 
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Figure 11. The cumulative mass within R (M T ) as a function 
of distance (R) from the center of a MSC for FG stars (red solid) 
and SG ones (blue dotted) in the standard model. These FG and 
SG stars are old and new ones, respectively, in the present study. 



1 


1 




| , I , I | 

s rol -0.0.p lh -0 


0.8 






s rol -°-°.Pt h -i° 4 : 








; 


0.6 






s rol -0.06,p lh -10* 


0.4 


- \ 










"\ 




0.2 

















5 10 15 



Radius (pc) 

Figure 12. The mass fraction of SG stars (Fsg = Msg/(Mfg + 
Msg)) as a function of distance from the center of a MSC (R) 
for different four models with s r ot = and p tri = atoms cm -3 
(black solid), s ro t = and p t h = 10 4 atoms cm -3 (blue dotted), 
•Srot = 0.32 and p tri = atoms cm -3 (cyan short-dashed), s ro t = 
0.32 and p tri = 10 4 atoms cm -3 (green long-dashed), and s ro t = 
0.08 and p tri = 10 4 atoms cm" 3 (red dot-dashed). The initial 
masses of MSCs in these models are the same as that of the 
standard model (M s = 1O 6 M ). 

10 5 Mq in this parameter study) above which AGB ejecta 
can be retained effectively and (ii) R TC t is larger for larger 
M s . The reason of these dependences is that MSCs with 
larger M s have deeper gravitational potentials so that they 
can retain AGB ejecta more effectively. Differences in F Tet 
between R < lOpc and R < 50pc are small for models with 
Ms ^ 6 x 10 5 M Q , which means that most the accumulated 
gas can form very compact gaseous regions within the cen- 
tral regions of MSCs. Thus MSCs with larger M s can obtain 
a larger amount of AGB ejecta for further star formation: 
M s is a key for determining whether MSCs can have SG 
stars. 



4.2 Secondary star formation 

Fig. 7 shows how gas ejected from AGB stars is accumulated 
within the central region of a MSC and consequently con- 
verted into new stars (i.e., SG ones) there within 13.6 Myr 
in the standard model with M s = 10 6 M Q , i? s = 44.7pc, 
Srot = 0.0, and p th = 10 4 atoms cm -3 . The gaseous ejecta 
from AGB stars initially located in the central region of 
the MSC can be accumulated there within ~ lMyr so that 
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Figure 13. The final distributions of old stars (magenta) and 
new ones (cyan) projected onto the x-y plane (left) and the x- 
z one (right) in the more massive model with M s = 10 7 Mq, 
-Ra = lOOpc, Srot = 0.02, and p t h = 10 4 atoms cm -3 . 
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Figure 14. The same as Fig. 11 but for the model with the more 
massive one shown in Fig. 13. 



star formation can start efficiently. The star formation from 
AGB ejecta can proceed like a "starburst" with most gas 
being consumed rapidly within the first 3 Myr evolution of 
the MSC. The star formation rate (SFR) reaches its maxi- 
mum value (O.O3M0 yr" 1 ) at T = TO Myr and then rapidly 
declines to be 0.001M Q yr" 1 ) at T = 6 Myr. 

About 86% of the gas can be converted into new stars 
within 13.6 Myr in this model. This very high star forma- 
tion efficiency (e s f ) results from the fact that almost all of the 
AGB ejecta can be accumulated very quickly (in less than a 
few 10 6 yr) within the central region of the MSC owing to 
the deep gravitational potential and consequently converted 
into new stars there. The accumulated gas is strongly self- 
gravitating (i.e., gas mass comparable to stellar mass there) 
in the central region of the cluster so that the gas can con- 
tinue to collapse to form new stars. This result clearly sug- 
gests that secondary star formation within MSCs is respon- 
sible for the origin of the observed high stellar densities of 
the present GCs: A GMC can not be converted into a GC 
just by one star-formation event. 

A new compact star cluster embedded in low-density 
residual gas can be formed in the central regions of the MSC 
within 13.6 Myr in this model. Owing to the high star for- 
mation efficiency (> 0.5), the new cluster is highly likely to 
survive from gas removal by energetic winds of massive stars 
and supernova explosion. As shown in Fig. 8, the new clus- 
ter is much more compact than the original MSC and has a 
half-mass radius (~ 3pc) significantly smaller than that of 
the original MSC (~ llpc). Owing to its compactness, the 
new cluster is much less susceptible to tidal destruction by 
its host galaxy. Thus a "nested stellar system" composed of 
a diffuse original cluster (FG stars) and a compact new one 
(SG ones) can be formed as a result of secondary star forma- 
tion from AGB ejecta within the MSC in this model. This 
nested structure is one of common features of the simulated 
MSCs in the present study. 

Fig. 9 shows how the fractions (F ns ) of the total mass of 
new stars (M ns ) to the initial total mass of the MSC (M s ) 
depend on M s for R < lOpc and 7? < 50pc in the models 
with p t h = 10 4 atoms cm -3 . Clearly there is a threshold 
cluster mass (M t h) above which secondary star formation is 
possible: M t h ~6x 10 5 M Q and does not depend on pth in 
the present study. Furthermore, F ns is higher in MSCs with 
larger M s , which implies that e s f is also higher in MSCs with 
larger M s . Almost all stars can be formed within half-mass 
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Figure 15. The radial profiles of rotational velocity V (solid) 
and one-dimensional velocity dispersion a (dotted) for FG stars 
(red) and SG ones (blue) in the more massive model shown in 
Fig. 13. Here R is the projected distance along the x-axis and the 
y-component of velocity for each star is used for estimating the 
line-of-sight velocity and velocity dispersion. These FG and SG 
stars correspond to old and new stars, respectively, in the present 
study. 



radii (~ lOpc) of MSCs with M s > Mth so that there can 
be no differences between F ns for R < lOpc and R < 50pc. 

Although e s f is rather high (0.86) in the model shown 
in Figs. 7 and 8, the total mass of new stars can be at 
most 4 x 1O 4 M0 owing to the adopted /agb = 0.04 that 
is reasonable for the standard IMF. This means that MSCs 
with M s ~ 1O 6 M0 can become low-mass GCs with multi- 
ple stellar populations. Furthermore, e B f is 0.3 for the model 
with M s = 6 x 1O 5 M0, which suggests that the new cluster 
can significantly expand after gas expulsion due to energetic 
winds of massive stars and supernova. The cluster is likely 
to become a low-mass, low-density GC with a smaller frac- 
tion of SG stars (like Pal 5). The dependence of e s f on M B 
implies that more massive GCs have larger fractions of SG 
stars. 

The details of star formation histories depend weakly 
on s ro t and pth such that gas can be more rapidly converted 
into new stars in MSCs with smaller s rot and p t h- Fig- 10 
shows that overall star formation rate is higher in the model 
with s r ot = (e S f = 0.86) than in that with s ro t = 0.08 (e s f = 
0.47). This is mainly because initial angular momentum of 
the MSC in the model with s ro t = 0.08 can suppress the 
formation of a very compact gaseous region in the center of 
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Figure 16. The same as Fig. 2 but for the model with the 
external gravitational field of the host galaxy of a MSC with M s = 
2 X 1O 5 M0. The MSC is initially located in the center of its host. 
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Figure 17. The same as Fig. 6 but for the models with external 
gravitational fields of hosts. The MSCs are initially located in the 
centers of their hosts. 



the MSC. The secondary peak of the star formation rate in 
the model with s rot = 0.08 at T = 2.8 Myr is due to later 
infall of gas onto the MSCs center where a compact gas disk 
is formed. 



4.3 Structures and kinematics 

The formation of nested structures with very compact new 
clusters and initially more diffuse old ones is one of essential 
ingredients of the present simulations. Fig. 11, which shows 
a typical example of the nested structures, describes how 
the total stellar mass within R depend on R for old and 
new stars in the standard model. Clearly the final cluster 
is dominated by new stars (SG) within R < 2pc, and the 
faction (Fsg) of the total mass of new stars (Msg) to the 
total mass of old stars (Mfg) and new ones is more than 
0.5 for R < 2pc. Fig. 12 shows how Fsg depends on R for 
different models with different p t h and s ro t- Irrespective of 
these model parameters, Fsg is higher in the inner regions 
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Figure 18. Star formation histories of low-density MSCs M B = 
1O 6 M0, Rs = 75pc, Srot = 0.0, and p t h = 10 4 atoms cm -3 in two 
models with (blue solid) and without (red dotted) external grav- 
itational fields of hosts. Normal and nuclear MSCs corresponds 
to models with and without external gravitational fields of hosts, 
respectively. 



of MSCs and lower in the outer ones. Fsg can be higher for 
higher p th and smaller s sot . 

As discussed in previous sections, M s should be at least 
2.5 x 10 6 M Q to explain GCs with M SG ~ 1O 5 M (for / AG b = 
0.04). Given that e s f is not 1 (0.3 — 0.9 for most models) 
and later dynamical evolution can reduce the total masses 
of MSCs, stellar structures and kinematics of the simulated 
MSCs with rather large M s (> 5 x 1O 6 M ) can correspond 
to those of the observed ones of GCs with masses larger than 
~ 10 Mq. Fig. 13 shows the final spatial distributions of old 
and new stars in the model with M s = 10 7 Mq, R b = lOOpc, 
Srot = 0.04, and p t h = 10 4 atoms cm -3 . Clearly the compact 
new cluster in the center of the MSC is significantly flattened 
owing to gaseous dissipation during gas accretion onto the 
central region. This significantly flattened shape can be seen 
in new stars for most models in which s rot is not (i.e., even 
the models with s ro t = 0.003 show such flattened shapes). 

The total mass of new stars (Msg) is larger than 1O 5 M0 
and the effective radius (Rh,SG) of the stars is 5.2pc, which is 
about a factor of 5 smaller than that of the old ones (J?h,PG)- 
The half-mass ratios (= Rh,sa/ Rii,fg) ar e very similar (~ 
0.2) in the present models with different model parameters. 
As shown in Fig. 14, the new stars dominate the inner region 
(R < 2pc) of MSC, though Fsg there can not become so 
high as that in the standard model. It is confirmed that 
more massive models with M s = 10 7 Mq with different s ro t 
and pth also show distinct nested structures. 

Fig. 15 shows clear differences in stellar kinematics be- 
tween old and new stars (thus FG and SG ones, respectively) 
in the final MSC for the model. New stars clearly have ro- 
tational kinematics with the maximum rotational velocity 
(Vrot) of 12.0 km s _1 and the smaller central velocity disper- 
sion (cto) of 9.5 km s _1 , which means Vrot/fo = 1.25. This 
rotational kinematics can not be seen in old stars, which 
have V r ot/o"o = 0.1 in this model. The models with larger 
Srot show larger Vrot/oo in new stars, which has already been 
discussed by Bekki (2010). These results suggest that the 
present GCs originate from MSCs composed of two stellar 
populations with significantly different kinematics. 
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4.4 Influences of host galaxies 

Fig. 16 shows the time evolution of the projected spatial 
distribution of AGB ejecta for a nuclear MSC in the model 
without star formation yet with the external gravitational 
field of its host galaxy dominated by dark matter halo 
(Mh = 10 9 Mq). In this model, only one gaseous sphere 
of a AGB star is included and located initially in the cen- 
ter of a nuclear MSC with M s = 2 x 10 5 M o (s ro t = 0.0) 
so that the results can be compared with the model with 
M s — 2 x 10 5 M Q shown in Fig. 2: It is much better to use 
these two comparative models with only one gaseous sphere 
of a AGB star, because the effects of MSC's host galaxy 
can be more clearly demonstrated. As shown in Fig. 16, al- 
though the gas sphere expands initially to a larger radius 
(R > 50pc), the gas can finally return back to the central 
region of the MSC within ~ 10 Myr owing to the deep po- 
tential well of its host galaxy. 

Fig. 17 shows that F rct in the models with external grav- 
itational fields of hosts do not depend so strongly on M B as 
those without the fields do: Irrespective of M B , most of AGB 
ejecta can be retained within the central 50pc of MSCs in the 
models with external gravitational fields of their host owing 
to their much deeper gravitational potential wells. However, 
F rct within R = lOpc depends more strongly on M B such 
that it is larger for larger A/ s . This results indicates that 
inner gas densities (R < lOpc) can not be so high in nuclear 
MSCs with lower M B so that star formation can not occur 
efficiently in the MSCs. It is actually confirmed that star 
formation can not occur efficiently in the low-mass models 
with M B s: 2 x 10 5 M Q , even if external gravitational fields 
are included. Thus the original masses of nuclear MSCs are 
still important for secondary star formation processes even 
if they are located in nuclear regions of their host galaxies. 

External gravitational fields can enhance star formation 
of nuclear MSCs, in particular, for those with larger M B 
and lower mass-densities. This enhancement can be most 
clearly seen in Fig. 18 which describes the time evolution 
of star formation rates for two models with M B — 1O 6 M0, 
R B — 75pc, Srot = 0.0, pth = atoms cm -3 , and with 
and without the external gravitational field from a host 
with Mh = 5 x 1O 9 M0. Star formation can start earlier 
and proceed more efficiently in the model with the exter- 
nal gravitational field (e SI = 0.78) than in the model with- 
out (e s f = 0.34). These differences between two models with 
and without external gravitational fields are less remarkable 
for high-mass and high-density MSCs, where star formation 
can proceed efficiently even in isolation. 



5 SURVIVAL OF ORIGINALLY MASSIVE 
STELLAR SYSTEMS 

In the present scenario, original MSCs need to keep sub- 
stantial masses at least for ~ 10 s yr after they form within 
GMCs so that gaseous ejecta from more massive AGB stars 
can be retained in MSCs and finally converted into new stars 
(=SG ones). The hosts (=massive dwarfs) can strip stars 
from MSCs owing to their tidal fields, even if the fields are 
not so strong as that of the Galaxy. Therefore we here inves- 
tigate how long MSCs can keep their substantial fractions 
of stars when they are under strong gravitational influences 
of their hosts. 
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Figure 19. Time evolution of the distribution of stars projected 
onto the x-y plane for models with M B = 1O 6 M0, R B = 75pc, 
fa = 1.2, M h = 5 X 10 9 M Q , r s = 1.6kpc, x s = 200pc, and 
f v = 0.75. Time T (in units of Gyr) that has elapsed since the 
simulation started is shown in the upper left corner of each frame. 
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Figure 20. Total stellar masses within lOOpc as a function of 
time T for models with r s = 1.6kpc and t„ = 1.0 (blue solid), 
r s = 1.6kpc and t a = 1.2 (green dotted), r B = 0.8kpc and t„ = 1.0 
(magenta short-dashed), r B = 0.8kpc and t a = 1.2 (red long- 
dashed). For these for models, M B = 1O 6 M M h = 5 X 1O 9 M 
and model parameters other than r B and t CT are fixed. 
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Figure 21. The same as Fig. 20 but for M B = 10 7 M C . 
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5.1 The model 

We consider how a MSC dynamically evolves after most of 
residual gas left behind from the formation of FG stars are 
expelled from the cluster through supernova events. Accord- 
ingly we adopt an assumption that (i) the MSC is assumed 
to consists purely of stellar particles and (ii) it is not neces- 
sarily in virial equilibrium owing to loss of the original mass 
of the GMC. The cluster with a mass M B and a size 7? s has a 
Plummer density profile with the scale length a s (R B = 5a s ). 
Owing to the above (ii), the central velocity dispersion a B of 
the MSC can be significantly larger than the one (<7 v ir) which 
the MSC can have when it is in virial equilibrium. Therefore 
we introduce a parameter f a that is the ratio of a B to cr v i r 
and ranges from 1.0 (virial equilibrium) to 1.3. Thus MSCs 
can expand significantly for models with f a > 1. 

The cluster is influenced by the host with the density 
distribution of the NFW halo (i.e., equation (7)). The initial 
location of the MSC with respect to the center of the host 
is set to be (x s , y B , z B ). The initial velocities of the MSC 
in the x-, y-, and z-directions are set to be (u B , v B , w B ). 
The orbital plane of the MSC is coincident with the x-y 
plane for all models in the present study (i.e., z B = and 
w B — 0). The initial direction of the velocity vector of the 
MSC is only a parameter for the orbit of the MSC owing 
to the spherically symmetric distributions of the host and 
the MSC. We therefore assume that the initial direction is 
parallel to the y-axis toward the positive y (i.e., v B > and 
m s = 0). We assume y B — 0, and thus x B and v s are free 
parameters that determine the orbit of the MSC. 

We consider that the orbit of a MSC is not necessarily 
circular within its host and therefore introduce a free pa- 
rameter / v that is the ratio of u B to the circular velocity v c i r 
at x B and ranges from 0.5 to 1.0 (i.e., circular orbit). We 
mainly investigate two representative cluster models with 
M s = 1O 6 M and R s = 75pc and with M B = 10 7 M Q and 
R s = lOOpc for different model parameters Mh, fa, x B , and 
f v . Although we run numerous models, we mainly describe 
the results of the models with Mh = 5 x 10 9 Mq and r B — 1.6 
kpc (c = 13.9). 

For each model, we count the total number of stars that 
are located within lOOpc from the center of the MSC at each 
time step in order to investigate the time evolution of the 
total stellar mass of the MSC. Stars that are once stripped 
and by accident located within the central lOOpc without 
being bound by the MSC can be counted as those in the 
MSC. We consider that since the number fraction of these 
stars is very small, these stars can only slightly overestimate 
the total stellar mass of a MSC. The initial stellar mass 
for a MSC is (1 — f a j)M B so that we can estimate the time 
evolution of the stellar mass for the MSC. 



5.2 Results 

Fig. 19 shows that the host of a MSC can gradually remove 
stars from the outer part of the MSC owing to the strong 
tidal field in the model with M B = 10 6 M Q , R B = 75pc, / CT = 
1.2, M h = 5 x 10 9 M Q , r B = 1.6kpc, x B = 200pc, and f v = 
0.75. About 60% of original stars can be stripped within 
~ 1 Gyr and the stripped stars are distributed widely within 
the original orbital plane and thus regarded as "field stars" . 
Although the mass and the size of the MSC can become 



much smaller than the original ones, the MSC is still strongly 
bound at T = 1.09 Gyr. The stripping process is gradual in 
this model: only 25% of the original mass of the MSC can 
be stripped within 0.1 Gyr. Thus, AGB ejecta might well be 
retained with the MSC to finally converted into new stars 
(=SG ones) in this model. 

Fig. 20 shows how the time evolution of total stellar 
masses within lOOpc of MSCs depends on f a and r B for 
models with M s = 10 6 M Q , Rs = 75pc, M h = 5 x 1O 9 M , 
x B — 200pc, and f v = 0.75. It is clear from this figure that 
stars can be more rapidly stripped from MSCs in the mod- 
els with larger f a (=1.2) and those in more compact hosts 
(r s = 0.8kpc). The MSC in the model with f a = 1.2 and 
r s = 0.8kpc can lose ~ 60% of its original mass within ~ 0.1 
Gyr owing to the combined effect of the more rapidly ex- 
panding stellar system and the stronger tidal field of the 
host. In this mode, AGB ejecta would not be so efficiently re- 
tained and consequently converted into new stars within the 
MSC. The mass fraction of stripped stars (Fstrip) is larger 
for models with smaller r B and larger f a . 

Fig. 21 shows how the time evolution of total stellar 
masses within lOOpc of MSCs depends on f a and r B for 
models with M B = 1O 7 M , R B = lOOpc, M h = 5 x 1O 9 M , 
x B — 200pc, and f v — 0.75. The results shown in Figs 20 and 
21 indicate that stars can be more slowly and less efficiently 
removed from MSC with larger M s owing to their stronger 
gravitational potential wells. MSCs with M B = 1O 7 M can 
not lose more than 50% of their original stellar masses within 
~ 1 Gyr so that AGB ejecta can be retained within the 
MSCs and converted into new stars. These results imply 
that original masses of the present GCs with multiple stel- 
lar populations should be rather large to keep substantial 
fractions of their original stellar masses to be immune from 
tidal destruction by their host within ~ 0.1 Gyr. 

Initial densities of MSCs and x B also can determine 
i^strip for a given set of parameters M B , f a , Mh, r B , and 
/ v . MSCs with lower densities and smaller x B are more sus- 
ceptible to tidal destruction of their hosts so that they can 
lose more stars within 1 Gyr (i.e., larger F at rip). Thus the 
time evolution of stellar masses of MSCs depends strongly 
on a number of parameters (e.g., M s and f a ) so that the 
time scale within which gaseous ejecta from AGB stars can 
be retained and converted into new stars can also depend 
on the parameters. The present parameter study, however, 
indicates that original MSCs are unlikely to be destroyed 
completely within a timescale of ~ 10 s yr for most mod- 
els: most MSCs can slowly (> 1 Gyr) disintegrate owing to 
the strong tidal fields of their hosts. We thus suggest that 
secondary star formation from AGB ejecta within MSCs is 
highly likely for most MSCs with M s = 10 6 - 10 7 M Q and 
can last at least ~ 10 s yr. 



6 DISCUSSION 

6.1 A threshold cluster mass for secondary star 
formation 

We have shown that M s can determine whether AGB ejecta 
from FG stars of MSCs can be retained within the MSCs 
and consequently converted into new stars efficiently. Fig. 
22 briefly summarizes the present key results on how M 3 
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Figure 22. Dependences of star formation efficiencies e s f (= 
M nB /M g , where M g is the initial total gas mass of a MSC) on M B 
for s ro t = and p th = 10 4 atoms cm' 3 . 



controles e a f for the secondary star formation. MSCs with 
larger M s can have higher e B f and we have confirmed that 
his dependence does not depend on s rot and p th . MSCs with 
more rotation in FG stars are more likely to show smaller e s f . 
Secondary star formation is possible if M B is larger than 6 x 
10 5 M Q , though e s f is not so high. As long as M s > 1O 6 M , 
e B f can be larger than 0.5 so that new compact clusters, 
which would not disintegrate even after gas expulsion, can 
form within old clusters. 

Thus our results have shown that there can be a thresh- 
old cluster mass (M th ~ [6 — 10] x 1O 5 M ) above which sec- 
ondary star formation with high e s f can occur within orig- 
inal clusters. This M t h can be as low as ~ 2 x 10 5 M Q , if 
the original clusters are ~ 10 times more compact than the 
scaling-rclations described by the equations (5) and (6) im- 
ply. However e B t is at most 0.1 even for such compact clusters 
with lower masses. On the other hand, even if MSCs with 
M s — 6 x 10 Mq are ~ 10 times more diffuse than the 
scaling-relations described by the equations (5) and (6) im- 
ply, their e s f can be as large as 0.1. Therefore, it is reasonable 
to claim that M th is (6 - 10) x 1O 5 M . 

D'Ercole et al. (2008) showed that secondary star for- 
mation within original clusters is possible for M B — 1O 7 M 
but not for M s = 10 6 M Q . This implies that M th can be 
somewhere between 1O 6 M0 and 10 7 Mq and is significantly 
larger than the derived Mth in the present study. This is 
mainly because modeling of star formation and treatment 
of AGB ejecta are different between these two works. Both 
models in D'Ercole et al. (2008) and the present study 
have advantages and disadvantages in describing evolution 
of AGB ejecta within clusters. For example, the present 
model enables us to investigate self-consistently 3D hydrody- 
namical interaction of numerous gaseous spheres from AGB 
stars whereas the model by D'Ercole investigated cooling of 
initially warm AGB ejecta self-consistently. More sophisti- 
cated models in our future papers need to investigate M th in 
a fully self-consistent manner by including physical effects 
that are treated by the present model in a simple way. 



6.2 Star-to-star abundance variation dependent 
on cluster masses 

We here discuss how M s determines nature of multiple stel- 
lar populations in GCs owing to existence of M t h- We pro- 
pose that M s is one of key parameters which can deter- 
mine differences in ages and various chemical abundances of 
light and heavy ones between multiple stellar populations. 
Based on the present numerical results, we divide MSCs into 
the following categories according to their M s . It should be 
stressed here that M B for a cluster is not the present mass 
(M gc ) but its original one. 



6.2.1 M B < 1O 4 M 

These low-mass clusters can neither retain AGB ejecta nor 
convert gas into new stars owing to their shallow gravita- 
tional potential wells. Therefore these clusters are highly 
unlikely to show abundance spread even in light elements. 
Observations showed that the Galactic open cluster (OCs) 
show unimodal distributions of CN band strength (Martell 
& Smith 2009) and do not have abundance spread in light el- 
ements (e.g., Carretta et al. 2007) and O-Na anti-correlation 
(De Silva et al. 2009). These low- mass clusters can be eas- 
ily destroyed by interaction with GMCs (e.g., Gieles et al. 
2006) so that they can not obtain fresh cold gas for further 
star formation. The observed low-mass OCs in the Galaxy 
and the Magellanic Clouds belong to this class. 



6.2.2 1O 4 M < Ms < 2 x 1O 5 M 

These clusters can not retain AGB ejecta efficiently within 
their central regions so that secondary star formation can 
not occur. However, they can obtain fresh cold gas for fur- 
ther star formation by capture and accretion of cold molec- 
ular gas that are either those left behind from previous star 
formation events or those located close to them (e.g., Bekki 
& Mackey 2009). Therefore, if secondary star formation oc- 
curs within these clusters, they can only show multiple stel- 
lar populations with age differences yet no/little abundance 
spread even in light elements. The intermediate-age clusters 
in the LMC and the SMC belong to this class. 

6.2.3 2 x 1O 5 M < Ms < 6 x 1O 5 M 

Clusters with this mass range can retain most of their AGB 
ejecta within the central regions but can not convert the 
stars so efficiently (e s f < 0.3) into new stars within a short 
time scale (< 10 8 ) yr. Although dynamical fate of the accu- 
mulated AGB ejecta for the cluster remains unclear, "trig- 
gering mechanisms" such as cluster mergers (or initially very 
high-density clusters) would be necessary to convert the gas 
so efficiently into new stars. Strong ram pressure stripping 
of the gas by ISM of their host galaxies could remove the 
gas from the clusters. 



6.2.4 6 x 1O 5 M < Ms < 1O 7 M 

The ordinary Galactic GCs with abundance spread in light 
elements can originate from clusters with this mass range, 
because AGB ejecta can be retained in the clusters and con- 
verted into new stars very efficiently (e s f > 0.5). The mass 
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fraction of SG stars in a present GC depends mainly on M s 
and dynamical evolution of its original cluster (i.e., destruc- 
tion by its host and long-term evolution driven by two-body 
relaxation) . If these massive clusters are initially stellar nu- 
clei of massive dwarf galaxies, then AGB ejecta from SG 
stars can be converted into new stars (i.e., the formation of 
third generation of stars). This sequential formation of stars 
in the nuclear regions of hosts can continue until the hosts 
are destroyed by much larger galaxies. 

6.2.5 10 7 M Q sj M s 

Very massive clusters with this mass range are progenitors 
of the present massive GCs with masses larger than 10 6 Mq. 
Stellar systems with M B > 10 7 M Q would be dwarf galax- 
ies themselves rather than clusters, and therefore at least 
some fraction of massive GCs (e.g., uj Cen and Gl) may well 
the remnants of nucleated dwarf galaxies. Owing to much 
deeper gravitational potential wells of nucleated dwarfs, not 
only AGB ejecta but also gas from more energetic massive 
stars and supernova would be able to be retained and finally 
converted into new stars. Thus such massive GCs can show 
abundance spread in both light and heavy elements. 

6.3 Mixing of AGB ejecta with fresh cold gas 

6.3.1 Bondi accretion 

Using analytical models, Pflamm-Altenburg & Kroupa 
(2009) recently have shown that if the mass of a cluster 
exceeds a threshold mass (= 10 6 Mq), then the cluster can 
obtain a significant amount of warm ISM from outside the 
cluster by Bondi accretion within a reasonable time scale 
(i.e., a few Gyr). We have confirmed their results using 3D 
hydrodynamical simulations, though the threshold mass for 
gas accretion is larger (= 3 x 10 Mq) in our estimation 
(Bekki et al. 2010). We suggest that although secondary for- 
mation from mixed gas of AGB ejecta and the accumulated 
ISM by Bondi accretion would be reasonable for massive 
GCs with possible age differences (like u> Cen), such star 
formation has the following potentially serious problems in 
explaining chemical abundances of ordinary Galactic GC. 

Firstly, if original clusters are within galaxies, then they 
can obtain gas initially located in different regions with pos- 
sible different chemical abundances (owing to radial and 
azimuthal abundance gradients with the galaxies) thus are 
highly likely to have gas and new stars with different [Fe/H]. 
Therefore, only a minor fraction of the Galactic GCs may 
well have SG stars formed from mixed gas of AGB ejecta and 
ISM obtained by Bondi accretion. Secondly, the Bondi accre- 
tion rate is too small (~ 10 _3 M Q yr -1 for M s = 3 x 1O 6 M ) 
for typical ISM so that an enough amount of SG stars can 
not form from gas chemically polluted only by massive AGB 
stars (> 5Mq) within an order of 10 8 yr. It is possible that 
if densities of warm ISM are as large as 10 2 atoms cm -3 , 
then clusters would be able to obtain an enough amount of 
gas within a timescale of less than 10 8 yr (but if this is the 
case, then such high-density gas may well be cold molecular 
clouds rather than warm ISM). 

Thirdly, even if warm ISM with initial temperature of 
~ 10 4 K can be accreted by clusters, further efficient radia- 
tive cooling (to 10-100K) within a short timescale (< 10 8 



yr) is required for further star formation. Our previous and 
present study have not yet investigated how warm ISM and 
AGB ejecta mixed with and consequently convert into cold 
molecular gas for further star formation using a sophisti- 
cated model, and accordingly we can not currently discuss 
the possibility of secondary star formation from the accumu- 
lated warm ISM. Our future more sophisticated models will 
investigate whether these three are really serious problems. 

6. 3. 2 Capture and accretion of cold molecular gas by 
clusters 

It is possible that the original MSCs would have a plenty 
of cold molecular gas left behind from the formation of FG 
stars, because the star formation efficiencies are highly likely 
to be well less than 100%. Given that GMCs which are pro- 
genitors of FG stars should be very massive (> 10 7 M Q ) in 
the present scenario, then the host GMCs would have nu- 
merous substructures (i.e., smaller molecular clouds) that 
remains intact even after the formation of FG stars. These 
molecular clouds can be captured by or accreted onto MSCs, 
as a previous simulation demonstrated already (Bekki & 
Mackey 2009). Owing to initially low temperatures of the 
captured cold gas, star formation from gas mixed between 
AGB ejecta and the captured gas would proceed efficiently. 

A potential problem of this process is that it remain 
unclear whether cold molecular clouds left behind the for- 
mation of FG stars can really keep the same abundances 
as those of FG stars (i.e., remain "pristine") without being 
chemically polluted by massive stars and supernova of FG 
stars. A significant fraction of the residual cold gas would 
be heated up to become high-temperature and low-density 
warm and hot gas so that they can not be captured later 
by MSCs owing to their high-temperature and low-density. 
Owing to clumpy structures of the original massive GMCs, 
a significant amount of energy from massive stars and su- 
pernova can be expelled from low-density inter-substructure 
regions of the GMCs toward the halos of their hosts. Some 
smaller cold molecular clouds that can survive from thermal 
and kinetic feedback effects of energetic stars would be able 
to be captured later by MSCs and consequently mixed with 
AGB ejecta. 

The above processes of mixing of cold molecular gas 
and AGB ejecta have not been investigated in detail by 
any previous 3D hydrodynamical simulations of star-forming 
GMCs. Bekki & Chiba (2007) investigated (i) how massive 
stars and type II supernova influence gas dynamics within 
star-forming GMCs and (ii) chemical abundances of gas left 
behind from star formation. They however did not include 
ejection of gas so that they could not investigate how AGB 
ejecta mixed with the residual gas to form new stars. It is 
thus doubtlessly worthwhile for our future more sophisti- 
cated models to investigate whether and what fraction of 
residual molecular cloud left behind from the formation of 
FG stars can be converted into new stars without being 
chemically polluted by massive FG stars. 

6.4 Origin of GCs with multiple stellar 

populations in the Magellanic Clouds 

A number of recent observational studies have found that a 
significant fraction of intermediate-age cluster in the LMC 
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and the SMC show double main sequence turnoffs (hereafter 
referred to as DMSTOs for convenience) on their color mag- 
nitude diagrams (Mackey & Broby Nielsen 2007; Mackey et 
al. 2008; Glatt et al. 2008; Goudfrooij et al. 2009; Milone 
ct al. 2009). The simplest explanation is that the observed 
DMSTOs in each of these clusters represent two distinct 
stellar populations with differences in age of ~ 100 — 300 
Myr. Recently Bastian & de Mink (2009) have proposed 
that stellar rotation in stars with masses between 1.2 and 
1.7M© can mimic the effects of DMSTOs without resorting 
to age differences between stellar populations within clusters 
of the Magellanic Clouds. However, the presence of a dual 
red clump of giant stars in the color magnitude diagram of 
NGC 419 in the SMC is suggested to be inconsistent with the 
claim by Bastian & de Mink (2009) but be explained nat- 
urally by the presence of multiple star-formation episodes 
(Rubele et al. 2010). 

Mucciarelli et al. (2008) investigated chemical abun- 
dances of light odd-Z, q, iron-peak, and neutron-capture 
elements for four intermediate-age clusters in the LMC and 
found negligible star-to-star scatter for them. However, Muc- 
ciarelli et al. (2009) revealed significant abundance inhomo- 
geneities in [Na/Fe], [Al/Fe], [O/Fe], and [Mg/Fe] and O-Na 
and Mg-Al anti-correlations in old GCs (> 10 Gyr). These 
observations imply that intermediate-age clusters have (at 
least) two generations of stars with different ages yet similar 
chemical abundances whereas old ones are like the Galactic 
GCs with abundance spread in light elements. Therefore, 
any theory of GC formation need to explain why the LMC 
appears to have two different types of GCs with multiple 
stellar populations. 

Bekki & Mackey (2009) showed that low-mass clusters 
(M s ~ 5 x 10 4 M Q ) can capture cold GMCs and thereby 
convert the gas into new stars within their central regions 
in the LMC. The present study has shown that such low- 
mass clusters can not efficiently retain AGB ejecta and thus 
are unable to use the ejecta for further star formation. The 
results by Bekki & Mackey (2009) and the present study 
therefore combine to suggest that if the observed DMSTOs 
in clusters of the LMC are due to age differences of stellar 
populations, then the two populations are highly unlikely to 
show abundance spread in light elements owing to their low 
masses (< 10 5 M©). 

The present study also suggests that the old clusters 
of the LMC can show abundance spread in light elements, 
mainly because their original clusters are significantly more 
massive than the present ones with masses ranging from 
1O 5 M to 6 x 1O 5 M . Recently Conroy & Spergel (2010) 
have suggested that clusters with masses larger than a 
threshold mass of - 10 4 M Q can form SG stars from AGB 
ejecta mixed with ISM accreted onto clusters. The observed 
lack of abundance inhomogeneity in light elements in the 
intermediate-age clusters of the LMC (Mucciarelli et al. 
2008) seems to be inconsistent with their model. The present 
study suggests that the threshold mass for the formation of 
SG stars from AGB ejecta mixed with pristine gas is much 
larger than the one suggested by Conroy & Spergel (2010). 

Old and metal-poor ([Fe/H]< —1.65) clusters in the 
LMC have masses more than 10° Mq (Mackey & Gilmorc 
2003) and would have lost a significant fraction of their 
masses due to long-term internal evolution effects and ex- 
ternal tidal field of the LMC. Therefore recent observations 



on star-to-star abundance variations in light elements in 
old and intermediate clusters of the LMC by Mucciarelli et 
al. (2008, 2009) strongly suggest that the threshold cluster 
mass (Mth) for secondary star formation from AGB ejecta 
mixed with pristine gas would be much more than 10 5 M Q . 
These observations are consistent with the present model 
with Mth ~ (6 - 10) x 10 5 M Q and the one by D'Ercole et 
al. (2008) which showed that clusters with M s ~ 1O 7 M can 
form SG stars from AGB ejecta. 

6.5 Were original massive stellar systems really 
"clusters" ? 

D'Ercole et al. (2008) first claimed that original single stel- 
lar systems ( "clusters" ) composed purely of FG stars should 
have masses of ~ 1O 7 M0 to explain the observed masses of 
massive GCs (~ 10 6 M Q ). Although our present simulations 
confirm their claim, it remain unclear how such massive sin- 
gle clusters form from very massive GMCs. Given the possi- 
ble substructures within host GMCs for FG stars, the orig- 
inal stellar systems can be still clusters of smaller clusters 
when FG AGB stars to eject their gaseous winds. If this is 
the case, evolution of AGB ejecta in such clusters of clusters 
would be significantly different from what D'Ercole et al. 
(2008) and the present study describe. 

It would be also possible that (i) original massive stellar 
systems of GCs are dwarf galaxies themselves and therefore 
(ii) compact stellar systems formed from AGB stars of the 
dwarfs are identical to their stellar galactic nuclei. In this 
scenario, stellar nuclei are dominated by SG stars that can 
be formed more quickly from AGB ejecta of FG ones ow- 
ing to much deeper gravitational potentials of their host 
with dark matter halos. Therefore, stellar nuclei of nucle- 
ated dwarfs are highly likely to be dominated by He-rich, 
Na-rich, and O-poor stars. Future observations on chemical 
abundances of stars in nuclei of the Sagittarius dwarf and 
NGC 205 will enable us to discuss whether this scenario is 
physically viable. 

6.6 Implications 

6.6.1 A bottom-heavy IMF for SG stars ? 

The present simulations have shown that secondary star for- 
mation from AGB ejecta can occur at the central regions of 
original clusters, where mean stellar number densities ex- 
ceed ~ 2 x 10 2 stars pc~ 3 (for M s = 1O 6 M ). Recently 
Krumholz et al. (2009) investigated formation of massive 
stars by gaseous accretion within a gas cloud with a mass 
of IOOMq and a size of 0.1 pc and found that two (binary) 
stars with masses of 29.2M© and 41. 5M© can be formed 
within 5.7 x 10 4 years. Therefore the present simulations 
imply that gas clouds for massive stars (as those investi- 
gated by Krumholz et al. 2009 above) can interact frequently 
with stars within the central regions of a dense stellar sys- 
tem owing to the small mean separation of stars (~ 0.2 
pc) and the short dynamical time scale (~ 5 x 10 5 yr). It 
would be accordingly possible that massive star formation 
is severely suppressed owing to this star-cloud interaction 
which can prevent gas accretion onto clouds forming mas- 
sive stars: IMFs for the formation of SG stars within MSCs 
would be bottom-heavy. If formation of binary star-forming 
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gas clouds, which are progenitors of binary stars, can be also 
prevented owing to star-cloud interaction, then the observed 
small binary fractions among SG stars in GCs (D'Orazi et 
al. 2010b) can be naturally explained. 

If only SG stars have a bottom-heavy IMF with the 
power-law slope of 3.35 with FG ones with a = 2.35, then 
the number fraction of low-mass stars with masses lower 
than O.5M for SG stars to that for FG ones is 1.7 for 
m u = 12OM0. This number fraction can become as large as 
2 by assuming unrealistically low m u (< IMq) and steeper 
a (< 3). Therefore, if IMFs can be bottom-heavy only for 
SG stars, then the required original masses of clusters that 
finally become the present GCs can be up to by a factor of 
2 smaller than those derived in the present work with the 
standard IMF both for FG and SG stars. Thus, the possible 
bottom-heavy IMF can not change the present main con- 
clusion that the progenitor clusters of the present GCs are 
much more massive than the present GCs. 

6.6.2 Origin of rotation, shapes, and scaling-relations of 
GCs 

Previous works discussed the origin of shapes and rotation 
(e.g., Frenk & Fall 1982; Einscl & Spurzem 1999), sizes 
and luminosities (e.g., Murray 2009), scaling-relations (e.g., 
Djorgovski 1993; Bekki et al. 2004; Gieles et al. 2010) using 
analytical models and collisionless numerical simulations. 
They accordingly did not discuss the origin of GC properties 
in the context of initially nested structures and ignored the 
importance of gas dynamics in the formation of GC prop- 
erties. If most GCs really originate from nested stellar sys- 
tems, then conclusions of the above previous works based 
on models with "non-nested" stellar structures must be dra- 
matically modified. Thus it would be reasonable to say that 
the above conclusions can apply only for low-mass clusters 
with single stellar population and no nested structures. 

Bekki (2010) recently has shown that the origin of the 
observed rotation in GCs (e.g., Meylan & Mayor 1986; Mey- 
lan & Heggie 1997; Anderson & King 2003) can be closely 
associated with dissipative gas dynamics of AGB ejecta from 
FG stars within forming GCs. Given the observed large mass 
fractions of SG stars in the present GCs (e.g., Carretta et al. 
2010), gas dynamics of AGB ejecta within original clusters 
composed only of FG stars may well be a key determinant for 
structure and kinematics of the present GCs. Although re- 
cent numerical simulations have investigated long-term dy- 
namical evolution of GCs with initially nested stellar struc- 
tures (e.g., Decressin et al. 2008; D'Ercole et al. 2008), their 
models do not include initial rotation and flattened shapes of 
SG stars so that their results can not be used for discussing 
the origin of rotation and shapes of GCs. Our future long- 
term (> 10 Gyr) dynamical simulations of GCs based on 
structures and kinematics of GCs with nested stellar struc- 
tures, flattened 3D distributions, and rotation simulated in 
the present study will enable us to discuss structure and 
kinematics in a consistent manner. 



6.6.3 Where are forming GCs with nested structures ? 

Recently Vinko et al. (2009) have found two stellar popula- 
tions with younger and older ones with ages of 10 — 16 Myr 



and 32 — 100 Myr, respectively, in the young, massive stellar 
cluster Sandage-96 in a spiral arm of NGC 2403. Given the 
observed possible large age difference (up to ~ 100 Myr), 
this type of young massive clusters can be the candidates 
that are now forming or have just formed SG stars from 
AGB ejecta from FG ones. Vinko et al. (2009) also have 
found that the younger population are located closer to the 
center of the cluster in comparison with the older one. If 
the cluster has a mass as large as 1O 6 M0, the observation 
would be consistent with the present results on the nested 
structure of MSCs, though quantitative comparison (i.e., the 
mass-ratio and the half-mass-radius-ratio of the two popula- 
tions) with the present results is not currently possible. The 
present study suggests that if the total mass of the cluster 
is as large as 1O 6 M0, the redder colors in the younger pop- 
ulation of the cluster can be due to gas and dust that now 
surrounds the younger population and were previous ejected 
from the older one. 

If future observations find young clusters with large 
masses (> 1O 6 M ) and two populations with the age dif- 
ferences as large as 10 8 yr in actively star-forming galax- 
ies, they can be forming GCs with nested structures and 
thus provide strong constraints on the formation models of 
GCs like the present one. As suggested above, new compact 
clusters can be shrouded by gas and dust left behind from 
the formation of SG stars until supernova events expel the 
gas and dust: the dust-shrouded new clusters would have 
very red color for their young ages and possibly they can 
not be even seen in optical bands if dust extinction is so 
heavy. If FG stars can not be clearly identified as compact 
clusters owing to their initially more diffuse distributions, 
then dust-shrouded new clusters within the FG stars might 
well be identified as giant HII regions with no optical coun- 
terparts. However it is not clear whether the observed gi- 
ant HII regions with no near-infrared cluster counterparts 
in starbursting luminous infrared galaxies with numerous 
super-star clusters (e.g., Alonso-Herrero et al. 2002) can be 
dust-obscured very young compact clusters (SG) embedded 
in older ones (FG). 

Massive stellar systems with masses of ~ 1O 7 M were 
discovered in ultra-luminous infrared-galaxies with dust- 
shrouded strong starbursts (e.g., Monreal-Ibero et al. 2007). 
Numerical simulations showed that such dusty starburst 
galaxies can evove into "E+A" galaxies with poststarburst 
populations with ages of 0.1 — 1 Gyr (Bekki et al. 2001). The 
present sutdy has shown that MSCs with M s = 1O 7 M0 and 
ages of less than 1 Gyr (before destruction of the old clus- 
ters) can have nested steller structurs. These results com- 
bine to imply that young GCs with nested structures are 
highly likely to exist in E+A galaxies. It is thus worthwhile 
for future observational studies to search for massive young 
clusters in E+A galaxies (e.g., NGC 5102) to provide a clue 
to the origin of multiple stellar populations of GCs. 



7 CONCLUSIONS 

We have performed 3D numerical simulations of star forma- 
tion from AGB ejecta in MSCs formed from GMCs within 
their host galactic building blocks. We have considered that 
stars initially in MSCs and those newly formed from AGB 
ejecta corresponds to FG and SG stars, respectively, in the 



Multiple stellar populations 19 



present GCs and thus discussed the origin of multiple stel- 
lar populations of GCs. The main results are summarized as 
follows. 

(1) Gaseous ejecta from AGB stars can be retained 
within MSCs with M s larger than ~2x 1O 5 M , if MSCs are 
isolated (i.e., not being influenced by ram pressure of ISM). 
The mass fraction of gas retained in a MSC to M s (F re t) 
depends on M s such that F TCt is larger in larger M s . The 
multiple hydrodynamical interaction of gaseous ejecta from 
numerous AGB stars can significantly increase gas mass ac- 
cumulated within the central regions of MSCs, which was 
first found by the present study. 

(2) If M s of MSCs exceed a threshold cluster mass M t h, 
gaseous ejecta of AGB stars can sink into the central regions 
and then form high-density gaseous regions so that new stars 
(i.e., SG stars) can be formed there with high star formation 
efficiencies (> 0.3). M th is demonstrated to be (6 — 10) x 
10 j Mq in the present study. Deep gravitational potentials of 
MSCs play a great role both in retaining the accreting AGB 
ejecta onto the central regions and in forming high-density 
gaseous regions there. Therefore secondary star formation 
within clusters is inevitable in MSCs with M s larger than 
M t h- Owing to the existence of M t h, young open clusters and 
super star clusters with M s < 6 x 1O 5 M can not evolve into 
stellar systems with abundance spread in light elements. 

(3) Owing to the high star formation efficiencies (e s f = 
0.3 — 0.9), very compact stellar systems can be formed from 
AGB ejecta within a time scale of ~ 10 7 yr (after massive 
AGB stars start to eject their winds). The half-mass radius 
of a new compact cluster is by a factor of 5 smaller than 
that of the old one for most models in the present study. 
The nested clusters (or "cluster within cluster") with more 
diffuse old cluster and very compact new one are character- 
istic for the present numerical models. 

(4) At most 1 — 4% of all stars can be from new stars 
formed from AGB ejecta for the canonical IMFs in the 
present study. Therefore, in order to explain the present to- 
tal mass of SG stars with masses of ~ IO^Mq in a GC, 
M s is required to be as large as (3 - 10) x 10 6 M Q . The 
required mass for the progenitor MSC for u> Cen can be 
as large as 10 8 Mq owing to its larger present mass. Given 
that e = 0.3 — 0.9 for secondary star formation, the original 
masses of progenitors (M B ) for the present GCs with masses 
of M gc need to be larger than 25M gc . 

(5) Structural and kinematical properties of new com- 
pact clusters formed within original old ones depend on s ro t 
such that new clusters can be more flattened and more 
strongly supported by rotation in MSCs with larger s ro t- 
The new clusters are much more compact and more strongly 
supported by rotation than the old ones in all models with 
initial rotation in old clusters. These structural and kine- 
matical differences between old and new clusters do not de- 
pend on model parameters such as M B , R B , and p th . Given 
that the present GCs are dominated by SG stars, the result 
implies that the physical origins of flattened shapes and in- 
ternal rotation of GCs are closely associated with formation 
processes of SG stars through gas dynamics within MSCs. 

(6) Evolution of MSCs initially in the very centers of 
their hosts can be significantly different from that of MSCs 
outside the centers, mainly because AGB ejecta can be more 
effectively retained within nuclear MSCs owing to much 
deeper gravitational potentials of their hosts. Irrespective 



of M s , most AGB ejecta can be retained in these nuclear 
MSCs. Star formation, however, can proceed within nuclear 
MSCs, only if M B exceeds a threshold value that is almost 
the same as M t h described above. Strong gravitational fields 
of hosts can enhance secondary star formation of nuclear 
MSCs. Multiple episodes of star formation (i.e., third and 
fourth generations of stars) are possible in nuclear MSCs 
with M s ^ Mth until their hosts are destroyed during accre- 
tion of the hosts on larger galaxies. The origin of multiple 
generations of stars observed in lo Cen and NGC 2808 can 
be closely associated with evolution of nuclear MSCs. 

(7) MSCs can lose gradually (~ 0.1 — 1 Gyr) significant 
fractions of their stars during dynamical evolution within 
their hosts owing to tidal stripping by their hosts. The strip- 
ping processes and the timescale depend strongly on M s , 
orbits, and initial locations of MSCs within hosts. If stars 
in a MSC enter into their AGB phases after the MSC has 
lost most of their stars because of tidal stripping, then the 
AGB ejecta can not be efficiently retained in the MSC ow- 
ing to its much shallower gravitational potential. Therefore 
gaseous ejecta from AGB stars with lower masses are less 
likely to contribute to secondary star formation owing to 
the destruction of the MSC by its host. Thus the destruc- 
tion of MSCs due to their hosts can provide a mechanism 
by which star formation can be truncated in MSCs. 

(8) Direct capture and accretion of cold molecular 
clouds by MSCs themselves can be much more efficient ways 
to dilute the AGB ejecta of MSCs (i.e., mixing of these 
gaseous components) in secondary star formation within 
MSCs. The molecular clouds need to have chemical abun- 
dances very similar to those of MSCs, and thus they are 
either (i) those that are located very close to MSCs or (ii) 
those left behind from the formation of FG stars. The ob- 
served hierarchical structures of GMCs are suggested to be 
responsible for the dilution processes during formation of 
SG stars. 

(9) Low-mass MSCs with M B < 10 5 M Q in the LMC can 
obtain cold molecular gas by capture of the gas but can not 
retain AGB ejecta owing to much shallower gravitational 
potentials of MSCs. Therefore the intermediate-age clusters 
with DMSTOs recently discovered in the Magellanic Clouds 
can have multiple stellar populations with age differences 
yet no abundance spread in light elements owing to their 
low masses. The initial masses of the old GCs in the LMC 
were significantly larger than their present masses of (1 — 
6) x 10 5 Mq so that they could retain AGB ejecta to form 
SG stars with chemical abundances different from those of 
FG ones. 

(10) It is suggested that formation of massive stars 
can be severely suppressed in very dense central regions 
of MSCs (i.e., a bottom-heavy IMF), because dynamical 
interaction between stars and forming molecular cores 
(where individual stars can form) can prevent the growth of 
the cores. If dense stellar systems are really environments 
within which only low-mass stars can form (owing to a 
bottom- heavy IMF), then the original masses of the present 
GCs can be to some extent lower than the proposed one 
(Ms > 25M gc ) in the present scenario. 

Thus GCs with the observed high central stellar den- 
sities (~ 10 4 stars pc~ 3 ) can not be formed directly from 
GMCs with single star formation events within GMCs in a 
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scenario explored in the present study. The origin of the ob- 
served high stellar densities of GCs can be associated with 
unique star formation processes within central regions of al- 
ready compact stellar systems (composed of FG stars). The 
deep gravitational potential wells of original clusters with 
much lower stellar densities (~ 10 2 stars pc~ 3 at half-mass 
radii) can play a role in enhancing dramatically star for- 
mation efficiencies and consequently forming new compact 
clusters. 

The progenitor stellar systems for the present GCs 
with multiple stellar populations need to be very massive 
(M s > 25M gc ) to retain effectively AGB ejecta from FG 
stars, dilute the gas with cold molecular gas, form new stars 
from the gas, and finally have the total stellar masses of SG 
stars similar to the observed ones in the Galactic GCs. The 
present scenario explains why open clusters in the Galaxy 
and intermediate-age GCs in the LMC do not show clear 
abundance spread in light elements. If GCs really originate 
from MSCs with large M s , then the host GMCs for such 
MSCs should be even more massive (10 7 -I0 8 M Q ). We plan 
to investigate how very massive GMCs can form and evolve 
into MSCs due to efficient star formation within massive 
dwarf galaxies at high z. 
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